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Abstract 



We consider here the effects of a non-vanishing cosmological term on the final fate 
of a spherical inhomogeneous collapsing dust cloud. It is shown that depending on the 
nature of the initial data from which the collapse evolves, and for a positive value of 
the cosmological constant, we can have a globally regular evolution where a bounce 
develops within the cloud. We characterize precisely the initial data causing such a 
bounce in terms of the initial density and velocity profiles for the collapsing cloud. In 
the cases otherwise, the result of collapse is either formation of a black hole or a naked 
singularity resulting as the end state of collapse. We also show here that a positive 
cosmological term can cover a part of the singularity spectrum which is visible in the 
corresponding dust collapse models for the same initial data. 



1 Introduction 



The gravitational collapse of a dust cloud has been studied extensively in the literature [T|], 
especially in the light of the cosmic censorship conjecture of Penrose 0. There are no 
globally regular solutions in these models, and a singularity always develops in an initially 
collapsing configuration. These singularities could be naked, or hidden behind an event 
horizon, depending on the nature of the initial data functions representing the density and 
velocity profiles for the collapsing cloud |I[ , from which the collapse evolves. It is also known 
that the shell-focusing singularities in dust models are strong curvature singularities ||. In 
this sense these are physically genuine singularities and the possibility of extension of the 
spacetime through the same does not arise. 

In investigations such as above, one generally uses the Einstein equations with a vanishing 
cosmological term. Recent observations however give us an evidence that a large part of the 
universe is possibly dominated by an energy component with negative pressure [§], § . Among 
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various possibilities available 0, one candidate that is often considered is a cosmological 
constant |7[|, A, or vacuum energy density corresponding to a positive sign of A. Such a 
cosmological constant would represent a spatially uniform energy density distribution, which 
is time- independent, and its positive value acts as a globally repulsive force field. Also, the 
cold dark matter models with a substantial component supplied by the cosmological term are 
among the models which best fit the observational data || . Therefore it is of value to revive 
the cosmological term, as a constant or even as a time varying quantity, in the Einstein 
equations. With such a perspective in mind we have investigated here the gravitational 
collapse, and structure of the singularity in the spherically symmetric dust models with the 
presence of a no n- vanishing cosmological constant. 

Dust models with a cosmological term are also known in the literature 0. These models 
can be matched with the Schwarzschild-de Sitter spacetime at the boundary of the cloud [|T0[] . 



But the studies so far have been restricted to special cases, whereas we analyze here the 



Lemaitre-Tolman-Bondi (LTB) collapse models 0] with a cosmological constant, which 
provide the general solution to Einstein equations with dust as the source term. While 
studying these models, the weak energy condition, i.e. T^V % V^ > 0, for all non-spacelike 
vectors V 1 , is assumed everywhere in spacetime for the matter, which would be a physically 
reasonable condition to assume for the case of a collapsing cloud. 

In the next Section 2 we introduce the basic model and analyze the Einstein equations 
to check whether we can have globally regular solutions, as shown to exist in the case of 



homogeneous density dust solutions with a positive cosmological constant |12j . We also derive 



the general solution to the Einstein equations in the case corresponding to the marginally 
bound case in the LTB models (i.e. the case when energy function / = 0), and discuss 
the condition for avoiding shell-crossings. In Section 3 we investigate the structure of the 
singularity by studying (O)outgoing (Ra)dial (N)ull (Ge)odesics (ORANGEs) near the same. 
The bearing of the initial conditions - in the presence of the cosmological term - on whether 
the final state results in a black hole or a naked singularity is then analyzed. Section 4 
provides the conclusions. 



2 Dust collapse with A term 



In this section we first study the basic set of Einstein equations and the regularity conditions 
for collapse. Then we will look for the possibility of regular solutions when there is a non- 
vanishing cosmological term present, wherein an initially collapsing cloud rebounces at some 
later epoch so that singularity does not form. This corresponds to the occurrence of three 
phases, which are, collapse, reversal and subsequent dispersal. Finally, we give here a general 
solution with non-zero A term for the marginally bound case. 

2.1 Basic equations and regularity 

The model for a self-gravitating, spherically symmetric, inhomogeneous dust cloud with 
cosmological constant is given by the metric, 

ds 2 = -dt 2 + R ' c/ - dr 2 + R 2 (d6 2 + sin 2 : 6d(f> 2 ), (1) 
l + /(r) 
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where (t,r,9,4>) is a comoving coordinate system. Here f(r) > — 1, and r and R are the 
shell-labeling comoving coordinate, and the physical area radius respectively. In our notation 
the dot and the prime denote partial derivatives with respect to t and r respectively. The 
equation of state for matter in the interior of the cloud is that of dust (corresponding to the 
approximation of a fluid with negligible pressures), and the stress-energy tensor is given by, 

where e(r, t) is the energy density of matter. We assume the matter to satisfy the weak energy 
condition, which implies that e(r, t) > 0; it is equivalent to the strong energy condition as 
the principal pressures are zero. The Einstein equations in the presence of a cosmological 
constant can be written as, 

* 2 = ^nT + m + 1 R2 ' (2) 

We will be mainly considering the collapse situation (i.e. R < 0). The dust cloud is 
characterized by two free functions in general, representing the total mass (F(r)), and the 
energy (/(r)), inside the shell labeled by comoving coordinate r. The cosmological term A 
can in principle be of either sign, however the recent observations as indicated above seem to 
favour the positive sign. The above equation for R can be in principle integrated, and after 
integration one gets a constant of integration. One can fix the constant of integration by 
using the scaling freedom. We fix this by setting R = r on the initial hypersurface (t — 0). 
The two free functions F(r) and f(r) can be fixed by prescribing the initial density and 
velocity profiles through (j^) and @ respectively. We assume these free functions to have 

5 



the form 

F(r) = Fqt 3 + F n r 3+n + higher order terms, 
f{ r ) = fo r2 + fm r2+ni + higher order terms, 

where the choice of the first non-vanishing term is made in order to have a regular initial data 
on the initial surface t = from which the collapse evolves |IJ . From (|3]) it is clear that we can 
have both shell-crossing (R'(t sc (r),r) = 0), and shell-focusing singularities (R(t s f(r),r) = 0) 
in these spacetimes, depending on the dynamics of the shells given by (fj). In this paper we 
will only consider cases where there are no shell-crossings, because these are generally not 
considered to be genuine singularities, and as our main interest is in studying the nature of 
physical singularity corresponding to R — where the matter shells shrink to zero radius. 
This puts some restrictions on our initial data which we will discuss later. 

The metric exterior to the collapsing cloud has the Schwarzschild-de Sitter form, 

ds 2 = -gdt 2 + g- l dr 2 + r 2 {dd 2 + sm 2 6d<f) 2 ), (5) 

where g = 1 — 2M/r — Ar 2 /3. Matching the solutions (0) and (|3p at the boundary r& of the 
collapsing dust cloud one obtains 2M = F(rb) Jl3|, [T^]. It is known that the boundary can be 



made to bounce from an initially collapsing phase by choosing appropriate initial mass and 
cosmological term, as we see bellow. This behaviour can also be understood by analyzing 
the curve of allowed motion (fj) (for details see, |L3|], fL2|). 
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2.2 Regular solutions and rebounce 



First we look for the possibility of having globally regular solutions due to the presence of 
a cosmological constant. The equation (g), governing the dynamics of collapsing shells, can 
be conveniently written as 

Here V(R, r), defined as 

V(R,r) = 3F(r) + 3f(r)R + AR 3 , (7) 

is an analogue of the Newtonian effective potential governing motion of the shells. The 
allowed region of motion corresponds to V > as R 2 has to be greater than zero. 

If we start from an initially collapsing state, we will have a rebounce if we get R = for 
a given shell before that shell becomes singular. This can happen only if the equation 

3RR 2 = V(R, r) = 

has two real positive roots. In what follows it is convenient to define a quantity ((t, r) = R/r. 
The equation (|7|) can be rewritten as 

V = 3F + 3f(r + AC 3 r 3 . (8) 

which is a cubic equation with three roots in general. /,From the theory of cubic equations, if 
all the three roots of the equation above are real then at least one of them has to be positive 
and at least one negative. Note that V(R = 0, r > 0) = 3F > 0. Hence, any regular region 
between R = and the first zero of eqn. (§), i.e., Ci > 0, always becomes singular during 
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collapse and so we need two real positive roots for the above equation (||) for a possibility of 
rebounce. The region between the two real positive roots is forbidden region. So if we start 
on the right side, the collapsing shells bounce back and then we have continuous expansion. 
Since one of the real root has to be negative and the region between two real positive roots 
is forbidden, it is not possible to have oscillating solutions in these spacetimes. 

In the case of A = it is well-known that we cannot have a rebounce and the collapse 
necessarily results in a singularity. The cubic then reduces to a linear equation and the 
solution is given as, R = R ma x( r ) = —F/f. So only in the case / < we can have R = for 
positive R. This corresponds to the maximum possible physical radius for a given shell i.e. 
even if we start with initial expansion a given shell will reach the maximum radius R m ax{ r ) 
and then it will recollapse. 

In the case when A < 0, one and only one root is always positive. The other two roots 
are negative if 9F 2 < — (4/ 3 /A), or else they are complex conjugates. Therefore any initial 
configuration becomes singular in this case. The real positive root in this case gives an 
upper bound on the radius R = R ma x( r ) of a shell labeled r. This upper bound occurs as the 
negative (attractive) contribution from the A term keeps on increasing with an increasing 
R, while the contribution from gravitational attraction keeps on decreasing and so at some 
point for any value / the attraction due to A starts dominating. Hence, even if we have an 
initially expanding configuration, finally we always must have collapse in this case. 

For A > and f(r) > one can easily see that we can never have ( 2 = R 2 /r 2 = 0, so 
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singularity always forms if we start from an initial collapse. 

For the case A > and /(r) < 0, one root of the cubic above is always negative. If 

F> > -§£, (9) 

then the other two roots are complex conjugates. So the singularity always forms in such a 
case in initially collapsing configuration. On the other hand, if the initial data is such that 

F 2 < 

9A' 

then the other two roots are real positive. Let us denote them by £i and ( 2 with (i < ( 2 . 
The region between the two roots is forbidden. The entire space of allowed dynamics is given 
by the two disjoint regions [0, Ci] and [£ 2 ,oo]. If the initial scale factor ( lies in the first 
section, i.e. if Co < Ci> then we always have the singularity as the end point of collapse. Here 
r(i represents the upper bound for the physical radius of a shell in this region. If ( lies 
in the second section, i.e. Co > C2, then we will have a rebounce from the initial collapsing 
configuration. After the rebounce the physical radius of the shell keeps on increasing forever. 
There is no upper limit for the maximum value of ( in this region and r( 2 gives the lower 
bound for the physical radius of a shell, i.e. the shell rebounces at R — r( 2 . 

^From the above discussion we can see that we can have rebounce only in the case when 
A > and / < 0, and when the following two conditions are satisfied, 

4f3 

< -4. do) 
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A 


/>o 


/ = o 


/<o 


A < 


Closed solution, 
No rebounce 


Closed solution, 
No rebounce 


Closed solution, 
No rebounce 


A = 


Open (hyperbolic) solution, 
No rebounce 


Marginal (parabolic) solution, 
No rebounce 


Closed (elliptic) solution, 
No rebounce 


A > 


Open solution, 
No rebounce 


Open solution, 
No rebounce 


Various cases occur 
as in table ^|. 



Table 1: The various types of dust solutions with A 



F 2 > 

r ^ 9A 



F 2 < 

r ^ 9A 



F 2 < 

r ^ 9A 



1 < 2( 



Co < Ci 



Y- + | arccos 



9F 2 A 
4/ 3 



Co > C2 
1 > -( - 1 f ' US 



1 / 9F 2 A 

g arccos w |jg- 



Open solution, 
No rebounce 



Closed solution, 
No rebounce 



Open Solution, 
Rebounce 



Table 2: The various types of dust solutions with A > and / < 
and if Co > C2- With our scaling, this later condition can be written as, 

/ 



1 > 2 



~ 2 A 



cos 



1 9F 2 A 
- arccos 



(11) 



The tables p] and |2| give a systematic classification of the various possibilities we discussed 
above. 

The physical quantities like central density, p c = F /^ 3 (t), and curvature scalars 

'C 4 + C 2 C 2 \ „„„ 1f1 /C 4 + CC 2 C + C 2 C 2 \ .„ ,/c 2 + cc\ 



R ijkl R ijk i = 12 I - '.r - I , R 13 R,a = 12 I ? : " " I , ft = 6 



c 4 



1 



c 2 



stay finite as ({t, r) > for the regular models. 



There are several features here which are worth noting and have interesting physical 
significance as far as the dynamics of collapse is concerned, and which illustrate the effects a 
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non- vanishing cosmological term may have towards determining the final fate of a collapsing 
cloud of matter. Firstly, with a negative value of the cosmological term all the solutions 
become closed and singularity always forms in a future even if we start with initial expansion. 
This is to be expected because such a value will only contribute in a positive manner to the 
overall gravitational attraction of matter and just acts as a constant positive energy field 
helping the collapse. Next, as illustrated by the first line of Table ^|, there is a range of 
initial data where the collapse necessarily ends in a singularity despite however large positive 
value of the cosmological term. This is contrary to the belief sometimes expressed that a 
positive cosmological constant can always cause a bounce provided it is sufficiently large in 
magnitude. Finally, it is clear from above that there can be a rebounce only if the initial 
density is sufficiently low for a given positive value of A. This is so because the cosmological 
term becomes dominant with increasing distances, and gravity dominates at higher densities. 
Thus as the cloud is more disperse (lesser densities but bigger size), more is the effect of the 
cosmological term. 

2.3 The f = case 

As shown in the previous consideration, while a bounce and regular solution occurs for a 
specific range of the initial data, for a majority of the regular initial data space the collapse 
results into a spacetime singularity where densities and curvatures blow up. While for A = 
case we know in detail the structure of this singularity, and when it will be naked or covered. 
We would like to understand here the effects of a non-zero A towards the structure of the 
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singularity forming in such a collapse. 



To study these features for the dust collapse models, we analyze now explicitly the case 
/(r) = 0. While we have chosen / = for simplicity and clarity of the considerations, one 
expects quite similar behaviour in other cases also. The equation @ can now be written as, 



(12) 



where t c (r) is an integration function which it represents the time at which a given shell, r, 
becomes singular, i.e. R(t s f(r), r) = 0. The positive or negative sign respectively corresponds 
to the expanding and collapsing branches of the solution. In what follows, we consider only 
the negative sign because we are considering clouds which are collapsing, with R < initially. 
The integral (O) can be written as an infinite series in R near the centre as below, 



t-t c (r) = — 



2 R 3 ' 2 



F{r) 



-l) m (2m- 1)!! / AR 3 



i 2 m (2m + l)m! \3F(r) 



(13) 



Using the scaling freedom in our solution, we set R(t = 0, r) = r which determines t c (r) as 



2 r 3 ' 2 



1+E 



-l) m (2m- 1)!! / Ar 3 



3 Mr) I ~i 2-(2m + l)m! \3F(r 



(14) 



^From the above expressions we get, 



F' 

R = ~-=;R 

3F 



rF'' 
3F 



1/2 



' 3F + AR 3 / 2 \ 
3F + Ar 3 / 2 J [r 



1/2 



(15) 



Therefore, the condition R' > 0, implying no shell-crossings, can be satisfied if F' > and 
1 — rF'/3F > 0. This means that the total mass inside a shell and the matter density are 
increasing and decreasing functions of r respectively, as we move away from the centre. The 
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weak energy condition, implying positivity of energy density, guarantees that mass is an 
increasing function of r. Also, for any realistic density distributions, it would be physically 
reasonable that the density is higher at the center, decreasing away from the center. Thus 
we shall work with decreasing density profiles. Therefore, no shell-crossings occur in the 
spacetime for our choice of initial data, before the occurrence of the central shell-focusing 
singularity. 

As we will show in the next Section, the behaviour of collapsing shells near the centre 
would depend only on the first non-vanishing derivatives of the density and velocity profiles 
near the same. Therefore, the local visibility conditions are unaffected by the boundary 
conditions such as the initial choices of the mass function and the actual value ofA. On 
the other hand, the global behaviour of the trajectories coming out from the singularity can 
change due to the addition of a cosmological term. 



3 The structure of the singularity 



In this Section we will analyze the nature of the central singularity. In what follows we use 



the scheme, as developed earlier by Joshi and Dwivedi []T5|| , which gives us a necessary and 
sufficient condition for the local visibility of the singularity. The main idea here is to see 
if we can have ORANGEs in the spacetime, meeting the singularity in their past with a 
well-defined real positive tangent vector in a suitable plane. 
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The equation for the ORANGEs in spacetime ([]]), for /(r) = is 



d ± = g. 

dr 



(16) 



For convenience this can be written in the (u, R) plane as, 



dR R' 



<F A 



du ar a 1 



R 3 



R? 



;i7) 



where u = r a , and a > 1 is a constant to be determined later. ^From equation (|T7|) it is clear 
that there are no ORANGEs from the non-central part of the singularity curve because the 
first term under square-root goes to — oo as t approaches t s f(r). This means that we cannot 
have ORANGEs as solutions to above equation near the non-central singularity, and hence 
it is only the central singularity which can be possibly visible. 



We define X = R/u, to check if we can have a well-defined tangent for the equation QT7D, 
at r = 0, R = in the limit of t approaching the singular epoch £ s (0). Using the l'Hospitals 
rule we get, 



Xn = lim — = lim = lim — ~—r \ 1 



! + ) = &T(X 0| 0) (18) 



where the subscript denotes the value of the quantities at u = 0. The constant a can be 



uniquely fixed by demanding that R'/r a 1 is non-zero finite |Tj|. In the case that we are 
discussing, R' /r a ~ l remains finite if we choose a = 1 + 2n/3 |15[, and the above equation 



I8f) can be written as, 



1 

a 



X 



nF n 1 



X 



0. 



(19) 



3F JX\ 

where A = F/u. If this equation has a real positive root Xq then one will have at least one 
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null geodesic coming out of the singularity at R = 0,u = with the root X = X as a 
tangent in (u, R) plane. 

When n < 3, we have a < 3 and therefore A = 0, and the above equation reduces to, 

x 3/2 = K (20) 



2F ov /l + AF /3 

which always has a real positive root. Apart from an additional A term, this equation is 
analogous to the corresponding equation obtained for the LTB models, and reduces to the 
same for A = 0. Therefore, when either the first or the second derivative of density is non- 
zero the singularity is always, at least locally, visible. The addition of cosmological term 
changes only the value of the tangent (X ) to the ORANGES, but not the visibility property 
itself. Thus the corresponding dust naked singularity spectrum is stable to the addition of 
a positive cosmological constant. 

In this case, like the earlier studies on LTB models, the smaller root will be along the 



apparent horizon direction and a family of geodesies will come out along this direction [[U| . 
In this case we put a = 3 and the first term in the equation (|T7|) blows up and second term 
goes to zero such that the product is F . 



For the value n = 3, which corresponds to the critical case in LTB models [17 , we get 
A = F . Introducing X = F x 2 and £a = F 3 /F^ 2 Jl + AF /3 the root equation becomes, 



2x 4 + x 3 - £ A :r + a = 0. (21) 

This equation is similar to the corresponding case in the LTB models, with the modification 
in the definition of £. From the theory of quartic equations, this equation admits a real 
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positive root for £a < icrit = — (26 + 15/\/3)/2. Therefore, for a given central density F , 
and the inhomogeneity parameter as given by F 3 , the naked dust singularity can be partly 
covered by a positive A, as we get an additional positive term in the denominator in £a- 
But it is interesting to note that, however large, a finite A term cannot completely cover the 
corresponding visible part of dust models. 

In a similar manner we can see that a negative A will open-up some covered part in the 
dust spectrum. (As such our calculation does not use the positivity of A anywhere, so they 
go through even for negative A.) 

For all values n > 3 we have Ao = oo, and we cannot have a real positive root in these 
cases and the final singularity is safely hidden behind the event horizon. 

4 Conclusion 

Studying gravitational collapse of a dust cloud with a non-zero value of the cosmological 
term we observed that black holes, naked singularities and even globally regular solutions 
can develop as the final outcome of collapse. Each of these outcomes is determined by the 
choice of initial parameters, given in terms of the density and velocity profiles of the cloud. 
Though for simplicity we restricted to f(r) = case while analyzing the structure of the 
singularity, the results can be extended to the general case. We have also shown here that 
it is possible to cover a part of the naked singularity spectrum, in the corresponding critical 
branch of solutions in the LTB models, with the introduction of a positive cosmological 
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constant. An interesting conclusion that emerges is that the existence of naked singularity 
remains stable to the introduction of a cosmological term in Einstein equations. These results 
are of interest in that they allow us to understand the implications of a non-zero A towards 
the final outcome of gravitational collapse, in view of the recent observational claims about 
a non- vanishing cosmological constant. 

We acknowledge the support from University of Basque Country Grants UPV122.310- 
EB150/98, UPV172.310-G02/99 and the Spanish Ministry Grant PB96-0250. 
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